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ABSTRACT 



Aims. Expanding HII regions and propagating shocks are common in the environment of young high-mass star-forming complexes. 
They can compress a pre-existing molecular cloud and trigger the formation of dense cores. We investigate whether these phenomena 
can explain the formation of high-mass protostars within an infrared dark cloud located at the position of G327. 3-0.6 in the Galactic 
plane, in between two large infrared bubbles and two HII regions. 

Methods. The region of G327. 3-0.6 was imaged at 450 yum with the CEA P-ArTeMiS bolometer array on the Atacama Pathfinder 
Experiment telescope in Chile. APEX/LABOCA and APEX-2A, and Spitzer/IRAC and MIPS archives data were used in this study. 
Results. Ten massive cores were detected in the P-ArTeMiS image, embedded within the infrared dark cloud seen in absorption at 
both 8 and 24 fim. Their luminosities and masses indicate that they form high-mass stars. The kinematical study of the region suggests 
that the infrared bubbles expand toward the infrared dark cloud. 

Conclusions. Under the influence of expanding bubbles, star formation occurs in the infrared dark areas at the border of HII regions 
and infrared bubbles. 
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1. Introduction 

A necessary prerequisite to star formation is the existence of 
dense cores of self-gravitating gas. How the dense cores form, 
evolve and give birth to new stars remains a matter of de- 
bate for high-mass star formation (e.g., Motte et al. 2007). 
Possible phenomena responsible for their formation include ex- 
ternal events such as expanding HII regions and propagating 
shocks, which are common in the environment of young high- 
mass star-forming complexes (e.g., Zavagno et al. 2007) . The 
ionised front of an HII region can compress a pre-existing 
molecular cloud and trigger the formation of dense cores (e.g., 
Purcell et al. 2009). Stellar winds or shocks from supernovae 
can also lead to the formation of similar objects (e.g., Koo 
et al. 2008). The Galactic Legacy Infrared Mid-Plane Survey 
Extraordinaire (GLIMPSE; Churchwell et al. 2006) detected 
many infrared bubbles in the Galactic plane close to active sites 
of high-mass star formation. The infrared emission tracing the 
bubbles was interpreted as a signature of PAH emission in the 
photo-dissociation regions (PDRs) at the borders of HII regions. 
The bubbles may also be inflated by strong stellar winds. 



* We thank the staff of the APEX telescope for their support. APEX 
is a collaboration between the Max Planck Institut fiir Radioastronomie, 
ESO and Onsala Space Observatory 



We present the first 450 fim map of massive protostar pro- 
genitors within the high-mass star-forming region near the hot 
molecular core G327. 3-0.6 (Wyrowski et al. 2006) at the pro- 
jected extremity of an infrared dark cloud (Fig. 1). At a distance 
of 2.9 kpc (Bergman 1992), G327. 3-0.6 is surrounded by two 
large scale infrared bubbles and two HII regions, as shown in 
Fig. 1. The objective of this letter is to (1) provide a census of 
the objects in an early phase of star formation that are embedded 
in the infrared dark areas and (2) investigate the relationship be- 
tween early phases of high-mass star formation and the influence 
of the expansion of nearby infrared bubbles and HII regions. 



2. P-ArTeMiS 450 //m continuum imaging 

2.1. Observations 

The region of G327. 3-0.6 was imaged at 450 //m in November 
2007 with the P-ArTeMiS bolometer array on the Atacama 
Pathfinder Experiment (APEX) telescope located at an altitude 
of 5 100 m at Llano de Chajnantor in Chile. The P-ArTeMiS de- 
tector is a 16 x 16-pixel prototype of the ArTeMiS large-format 
bolometer camera being built by CEA Saclay for APEX (e.g., 
Talvard et al. 2008). Three individual maps, corresponding to a 
total effective integration time of 1 .4 hours, were obtained with 
P-ArTeMiS at 450 fim toward the G327. 3-0.6 region using a 
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Right Ascension (J2000) 

Fig. 1. Overview of G327. Left: Spitzer/IRAC image at 8 /mi of the large scale infrared bubble with G327. 3-0.6 and RCW97 at the 
border. Right: P-ArTeMiS 450-/im emission. The contours are the 450-//m emission at levels of 4, 10, 20, 40, 50, and 90% of the 
140 Jy beam -1 peak flux. The rms is about 0.5 Jy beam -1 on the map. 



total-power, on-the-fiy scanning mode. Each of these maps con- 
sisted of a series of scans in either Azimuth or Right Ascension 
taken with a scanning speed of 16" s _1 . The cross-scan step be- 
tween consecutive scans was 2". 

The atmospheric opacity at zenith was monitored by taking 
skydips with P-ArTeMiS and was found to be between 0.6 and 
0.8 at A = 450 /mi, corresponding to precipitable water vapour 
amount between 0.5 mm and 0.7 mm. No dedicated pointing 
model was derived for P-ArTeMiS, but the pointing corrections 
proved to be very similar i.e., to within ~10" of the corrections 
appropriate for the APEX-2A instrument. Pointing, focus, and 
calibration measurements were achieved by taking both short 
"spiral" scans and longer on-the-fiy beam maps of Mars and 
Saturn. We estimated the absolute pointing accuracy to be ~5". 
The flux amplitude was to within ~ 30% uncertainty. As esti- 
mated from Mars maps, whose angular diameter was ~14", the 
main beam had a full width at half maximum of ~ 10" and 
contained ~ 50% of the power. Offline data reduction, includ- 
ing baseline subtraction, removal of correlated skynoise and 1/f 
noise, and subtraction of uncorrected 1/f noise using a method 
exploiting the high level of redundancy in the data, were per- 
formed with in-house IDL routines (Andre et al. 2008). 

2.2. Results and analysis 

The 450 /im emission originates in a large infrared dark cloud 
(2x1 pc 2 ) near G327.3-0.6 and at the south-eastern border of the 
RCW97 HII region (Fig. 1). Spitzer IRAC and MIPS images as 
well as LABOCA/APEX and APEX-2A archives observations 
were used in addition to the P-ArTeMiS data set. 

Ten submillimetre sources (SMM) were identified using 
a version of the Gauss-clumps program and the MRE multi- 
resolution program based on wavelet transforms (Starck & 
Murtagh 2006), which has been customised for continuum im- 
ages (see Appendix B in Motte et al. 2007). This method allowed 
us to first detect and then characterise SMM2, 4, 7, and 8, which 
were not identified by eye in the original 450 /im map (Fig. 1). 
We used MRE to complete a series of filtering operations and 
provided views of the image on different spatial scales. A cut- 
off angular size of 24" (i.e., 0.3 pc at 2.9 kpc), on which scale 
SMM2, 4, 7, and 8 were resolved, was estimated. All wavelet 
planes of the original image up to this scale were summed to 
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Fig. 2. Left: The image and contours represents the filtered P- 
ArTeMiS 450-//m emission where the emission over an area 
larger than 24" has been filtered. The contours are the 450-/mi 
emission at levels of 2, 4, 6, 8, 12, 16, 40, and 75% of the 101 
Jy beam -1 peak flux. The rms is about 0.4 Jy beam -1 on the 
map. Right: Zoomed view of the IRAC 8-/mi emission image. 
Contours represent the filtered P-ArTeMiS 450-/mi emission. 



create images such as those shown in Fig. 2. We note that the 
original image can be described by the sum of all wavelet views 
plus the smoothed image (last plane containing all remaining 
scales). Table 1 presents the SMM source characteristics {F pea k, 
the peak flux in the original map; F24" , the peak flux in the fil- 
tered map; FWHM, the full width half maximum size; Fj„ h the 
integrated flux) from both the filtered map (Fig. 2) and the orig- 
inal map (Fig. 1). The brightest source, SMM1, coincides with 
the position of the hot core and an infrared source (IRS1 in Fig. 
2). SMM4 coincides with a 8 and 24 /mi infrared source (IRS2 
in Fig. 2). Other sources do not have any infrared source coun- 
terpart in the MIPS or IRAC image. 

Assuming optically thin dust continuum emission at 450 /mi, 
the gas masses M of the sources were derived from the measured 
flux densities F450 Mm (either F or F24") using 



M = 



F45O umD 



*"450 urn #450 um(Jd) 



(1) 



where D is the distance to the source, K450 pm is the dust opacity 
per unit (gas + dust) mass column density at A — 450 /mi, and 
fi450 J um(r c i) is the Planck function B v {Ta) for the dust tempera- 
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ture 7d. We adopted K45o^ m = 0.04 cm 2 g _1 , which is consis- 
tent with the dust opacity law estimated by Hill et al. (2007) 
and appropriate in regions of moderately high gas densities ( 
«h, ~ 10 5 cm 4 ). The mean dust temperatures of SMM1, 2, 3, 
5, 6, 9, 10 were estimated from the observed 450 and 870 /mi 
flux density ratio by assuming optically thin dust emission with 
an emissivity index /3 — 2. The 870 //m flux density was derived 
from LABOC A/APEX observations. For SMM1, the hot core, 
and SMM 3, 5, and 9 in the RCW97 PDR, the temperature es- 
timates are probably lower limits, corresponding to upper mass 
limits. For three sources, SMM4, 7, and 8, the flux density ratio 
map did not provide any accurate measurement. A mean temper- 
ature for this part of the dark cloud was estimated to be ~ 20 K 
and used for these three sources (see Table 1). 

Mass estimates of the SMM sources were derived from both 
the peak fluxes and the integrated fluxes over their deconvolved 
sizes. There was an uncertainty of a factor of 2 on either side, 
mainly because of uncertain value of K450 SMM sources have 
masses that vary from 25 to 3800 M (Table 1). A second set of 
mass values, Mf™ AU , was also derived, which correspond to the 
estimate of the gas envelope masses of possibly dominant proto- 
stellar objects embedded within each SMM source. The diame- 
ter of the protostellar envelope is set to be 6000 AU (the typical 
fragmentation length-scale observed in high-mass star-forming 
regions - e.g., Longmore et al. 2006). The value of Mf^ AU was 
obtained by first calculating the mass corresponding to the peak 
flux density in the 10 beam and then applying a scaling factor 
to that mass, assuming a pocr~ 2 density distribution so that the 
measured mass scales as the aperture size. The value of M^ 0AU 
varies from 4 to 190 M (see Table 1). The overall mass of the 
dark cloud was estimated to be ~ 1 .7 x 10 5 M assuming a mean 
dust temperature of 20 K and taking into account the extended 
emission. 

SMM1 and SMM4 are associated with infrared sources at 
both 8 and 24 /mi (Fig. 2). Based on fits to approximate their 
spectral energy distribution from 24 to 870 /mi, the luminosities 
of SMM1 and SMM4 are estimated to be within the ranges of 
5xl0 4 - 10 5 L G and 3xl0 3 - 10 4 L , respectively. The luminosi- 
ties of the other SMM sources were not calculated because of 
the lack of mid-infrared emission at 24 /mi in the Spitzer/MIPS 
image. However, their luminosities were estimated to be within 
3xl0 2 - 10 4 L using their flux at 450 //m to extrapolate the 
luminosity values from those estimated for SMM1 and SMM4. 
The positions of SMM1 and SMM4 in an M em , - L bol diagram 
suggests that they are border-line Class 0/Class I objects that will 
evolve into massive stars of masses M» > 50 M and M* > 15 
M , respectively (Fig. 3). The other SMM sources, which have 
no detected IR emission, might represent earlier phases such as 
highly embedded class protostars or massive prestellar cores. 
Their mass estimates suggest that they contain enough material 
to form new high-mass stars. In particular, SMM3 might repre- 
sent a candidate massive prestellar object that has formed at the 
border of the RCW97 HII region. 

3. Triggered star formation in the dark cloud ? 

The MIPS 24 /mi emission image clearly exhibits two infrared 
bubbles (Fig. 4). The eastern one corresponds to the large bubble 
seen in the IRAC image (Fig. 1). It is centred on a OH/IR type 
star (IRAS 15502-5424) that is associated with circumstellar out- 
flows of blueshifted OH maser velocities of between -20 and -47 
km s" 1 (Te Lintel Hekkert et al. 1991). The star LSR velocity is 
estimated to be -8 km s . Interestingly, many molecular clouds 
and high-mass star-forming regions (SF075; RCW98) are ob- 
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Fig. 3. Envelope mass versus bolometric luminosity diagram 
comparing the locations of the protostellar sources of G327.6- 
0.3 (stars with error bars) with the positions of low-mass Class 
I (filled circles), Class objects (open circles), and high-mass 
protostars in NGC3576 (squares - Andre et al. 2008), as well as 
intermediate-mass protostars in NGC 2264 (triangles - Maury et 
al. 2009). Model protostellar evolutionary tracks, computed for 
various final stellar masses assuming the accretion/luminosity 
history described in Andre et al. (2008), are superimposed. 
Evolution proceeds from the upper left to the lower right as indi- 
cated by two arrows on each track, plotted when 50% and 90% of 
the final stellar mass has been accreted, respectively. The straight 
lines show two M em , - Lj, i relations marking the conceptual bor- 
der zone between the Class (M em , > M*/e) and the Class I 
(M eml < M*/e) stage (where e is the local star formation ef- 
ficiency). The dashed line is such that M em , oc L fco/ , while the 
dotted relation follows M env oc L9- 6 , as suggested by the accretion 
scenario adopted in the tracks. 



served around the bubble with LSR velocities close to the sys- 
temic velocity of G327. 3-0.6 (Russeil & Castets 2004; Urquhart 
et al. 2006). Two HII regions are also present near G327.3-0.6. 
To the north west of the infrared dark cloud is RCW97. To the 
south west, lies the HII region known as BBW29401, which 
is the brightest object in the Ha image at -48 km s (Russeil 
2003). The second largest 24 /mi infrared bubble encompasses 
BBW29401 (Fig. 4). 

The region of G327. 3-0.6 was imaged in 12 CO(3 - 2) us- 
ing the APEX-2A receiver on the APEX telescope (Wyrowski 
et al. 2006) . Using these archives observations, a first moment 
map and line spectra were created to determine whether the in- 
frared bubbles expand, compress the infrared dark cloud, and 
trigger star formation in G327.3-0.6. The velocity field around 
G327.3-0.6 can be divided, to a first approximation, across an 
east-west axis in two velocity colours: a blueshifted velocity 
(~ -50 km s _1 ) domain to the west, mainly located west of 
RCW97 (position D in Fig. 4); and a redshifted velocity (~ -46 
km s" 1 ) domain to the east (position C in Fig. 4). Along the 
eastern border of the infrared dark cloud traced by the curved in- 
frared filament, the velocity remains at the value of -46 km s _1 . 
The CO line toward the RCW97 HII region exhibits a double 
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Table 1. Dense cores detected in the G327 region. 



Fragment 


RA 


Dec 


r peak _ 




FWHM 


Pint 




AjFWHMd 
in! h 


iV1 env 


name 


J2000 


J2000 


[Jy/beam] 


[Jy/beam] 


[" x "] 


[Jy] 


[K] 


[M ] 


[Mo] 


SMM1 


15:53:08.6 


-54:37:09 


140 


101 


19 x 16 


440 


20 


3769 


184 


SMM2'' 


15:53:10.1 


-54:37:01 


- 


19 


16x 11 


35'' 


25 


197* 


23 


SMM3 


15:53:04.0 


-54:35:34 


24 


13 


25 x22 


132 


15 


2129 


45 


SMM4* 


15:53:11.5 


-54:36:46 




11 


20 x 11 


25'' 


20 


214* 


20 


SMM5 


15:53:01.4 


-54:35:20 
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" using the original map in Fig. 1 . 

* using the filtered map in Fig. 2. SMM2 needs confirmation. 
T d value is uncertain within 5 K. 

d with a factor 2 of uncertainty on either side of the mass value. 
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Fig. 4. Left: First moment map of the CO(3-2) line emission (colour image) superimposed on the MIPS 24-ytzm emission image 
(black and white image). Blue to red colours correspond to blue- and red-shifted velocities. The grey contours represent the P- 
ArTeMiS 450-ji/m emission. Right: Spectra taken at 4 positions in the first moment map: A corresponding to SMM4 in the infrared 
dark cloud, B to RCW97, C to the border of the infrared bubble, and D to west of RCW97 between SMM3, 5 and 9. 



profile with peaks at -50 and -46 km s (position B in Fig. 
4). Moving to the west and south-west of RCW97, the veloc- 
ity shifts toward -50 km s _1 . The morphology of RCW97 in the 
infrared also divides into two regions, with SMM3, 5, and 9 in 
between (Fig. 2). The CO spectrum at the position of SMM4 
exhibits multiple velocity components of velocity between -50 
and -40 km s _I (position A in Fig. 3). The 4 km s _I range of 
velocity peaks might indicate that the bubbles expand if they are 
not located at the distance of G327.3-0.6, but in its foreground 
and background. Similar profiles of CO spectra are modelled by 
Hennebelle et al. (2006) to explain a rapid increase in external 
pressure driving a compression wave in the Coalsack region, but 
of a far smaller velocity range. In any case, the CO line emission 
is probably self-absorbed by the cloud, which to some extent 
could also explain the spectra. 



The kinematical properties and the detection of the infrared 
bubbles at both 8 fim, a signature of PAH, and at 24 /mi, a signa- 
ture of heated interstellar dust, suggest that they possibly trace 
the interface between expanding fronts traced by PAH emission 
and the molecular cloud traced by thermal dust emission and gas. 



4. Conclusions 

The P-ArTeMiS image has detected massive cores embedded 
within the infrared dark cloud seen in absorption at both 8 and 
24 fim. Three SMM sources are detected at the western border of 
the RCW97 HII region. These SMM sources also coincide with 
infrared dark areas. The luminosity estimates of SMM sources 
indicate that the massive cores host luminous objects and there- 
fore form high-mass stars. 

The RCW97 HII region and the infrared dark cloud emit CO 
lines within the same velocity range. A possible interpretation is 
that both the infrared dark cloud and RCW97 are embedded in 
the same parent molecular cloud. Under the influence of expand- 
ing bubbles, star formation now occurs on the edge of RCW97, 
in the infrared dark cloud and at the border of HII region itself. 
In conclusion, circumstantial evidence of triggered star forma- 
tion has been observed toward G327.3-0.6. 
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